COMPARISON OF GEANT4 WITH EGSnrc FOR SIMULATION 
OF GAMMA-RADIATION DETECTORS BASED ON 
SEMI-INSULATING MATERIALS 

A.I. Skrypnyk*, A. A. Zakharchenko, M.A. Khazhmuradov 

National Science Center "Kharkov Institute of Physics and Technology", 61108, Kharkov, Ukraine 

(Received August 1 , 2011) 



O 

<N 
O 

Q 



We considered GEANTA version 4.9.4 with different Electromagnetic Physics Package for calculation of response 
functions of detectors based on semi-insulating materials. Computer simulations with GEANTA packages were run 
in order to determine the energy deposition of gamma-quanta in detectors of specified composition (Hgi2 and TIBr) 
at various energies from 0.026 to 3 MeV. The uncertainty in these predictions is estimated by comparison of their 
results with EGSnrc simulations. A general good agreement is found for EGSnrc and GEANTA with Penelope 
2008 model of LowEnergy Electromagnetic package. 

PACS: 29.40.Wk, 85.30De 
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1. INTRODUCTION 

The variety of electrophysical characteristics (spe- 
cific resistance, product of mobility /i and mean 
drift time r for holes and electrons - (fir)h, e ) is 
a serious disadvantage of semi-insulating materials 
(wide band-gap semiconductor with high resistivity) . 
Now this is a main cause that the semi-insulator 
gamma-radiation detectors could not be mass pro- 
duced. Wide band-gap semiconductor detectors have 
considerable spread of (/jLr)h,e values (in order of mag- 
nitude and more) [1] even if they are produced from 
one ingot. As result the same size detectors biased 
to the same voltage U & have a different charge collec- 
tion efficiency (CCE). It results in non-uniformity of 
serial devices response and in necessity of individual 
setup parameters selection. 

Also the material non-uniformities are a serious 
obstacle to comparison of detector characteristics 
measured in different studies. Computer simulation 
is an optimal method to overcome the material non- 
uniformity problem. Simulation may be useful both 
to research of wide band-gap detector behavior and 
designing of device based on them [2] . 

Previously we have developed the model of 
the planar wide band-gap semiconducting gamma- 
radiation detectors where EGSnrc Monte-Carlo (MC) 
simulation package was used for determination of the 
energy deposition of gamma-quanta and charged par- 
ticles [3] . This model was tested on the several groups 
of CdTe and CdZnTe detectors [4]. The observed dif- 
ference between model and experimentally measured 
amplitude distributions of radiation sources 137 Cs 
and 152 Eu was explained in process of this model ver- 
ification [2]. However, by this now all factors defining 
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such important gamma-radiation detector character- 
istics as sensitivity 5 and charge collection efficiency 
CCE are not determined even for the most researched 
semi-insulating materials as CdTe and CdZnTe. By- 
turn it does not permit to define exhaustive set of 
the control parameters of the wide band-gap detec- 
tor model. The determination of these parameters 
would allow calculating the most realistic response of 
the gamma-radiation detector. 

The model [4] enabled to obtain a good agree- 
ment between calculated and experimental response 
functions of CdTe (CdZnTe) gamma-radiation de- 
tectors in the most cases analyzed. However EGSnrc 
package possibilities are restricted only by a simula- 
tion of photon, electron and positron transport. It 
does not allow using the above mentioned model [4] 
for analysis of the charge collection efficiency experi- 
ments with the wide band-gap detectors when proton 
and a-particle beams are used [5]. Also the model 
[4] is unsuitable for researching of actual problems of 
radiation resistance of semi-insulators in mixed radi- 
ation fields (charged particles and/or neutrons and 
gamma-quanta) [6]. 

GeantA [7] is an actively developing toolkit for the 
simulation of the passage of charged particles, neu- 
trons and gamma-quanta through matter. This MC 
package offers a set of physical process models to de- 
scribe the interaction of charged and neutral particles 
with matter in the wide energy range. GeantA pro- 
vides various models of the same electromagnetic pro- 
cesses (EM packages) [8] . To implement the model [4] 
on GeantA platform properly it is necessary to define 
EM package which provides better agreement with 
EGSnrc simulation results. 

In the present work statistical characteristics of 
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the ^G^nrc-calculated response functions of planar 
wide band-gap detectors based on Hgl2 and TIBr 
to gamma-quanta with energies between 0.026 and 
3 MeV is compared with the results of GeantA simu- 
lation. The same radiation source and detector geom- 
etry parameters are used to calculate the photon and 
charged particle transport in EGSnrc and GeantA 
simulations. A good agreement between results of 
EGSnrc and GeantA simulations is obtained with us- 
ing Penelope-2008 EM package GeantA v.4.9.4 [9]. 

2. THE GEANTA AND EGSnrc 
ELECTROMAGNETIC MODELS 

G cant A includes three EM packages, which use differ- 
ent models for cross-sections of photons and charged 
particles interactions with the matter and different 
models for final state sampling algorithms. These 
are Standard EM model and two models for low en- 
ergy region identified as Livermore and Penelope 
EM model. All models describe the photoelectric ef- 
fect, the Compton scattering and the gamma con- 
version. The model processes of the induced emis- 
sion and scattering e~-e~ and e~-e + are available 
for electrons and positrons. 

In the recent years many articles were published 
where GeantA simulation results are compared with 
yield of preferences MC codes EGSnrc and MCNP 
[7-10]. EGSnrc and MCNP packages have longer 
application history in radiation physics and nuclear 
medicine in comparison with Geant4. As a result of 
often modifications of GeantA code and fixation of al- 
gorithm bugs the comparison of different GeantA ver- 
sion simulations with other Monte-Carlo codes and 
experimental data can show the differing results. In 
the present work GeantA v. 4. 9. 3 and v.4.9.4 simula- 
tions are compared with results of EGSnrc version 
4r2.3.1. 

The Standard model uses an analytical approach 
to describe the electromagnetic interactions in the 
range from 1 keV up to about 100 TeV. An analytical 
approach combines numerical databases with analyt- 
ical cross-section models assuming quasi-free atomic 
electrons while the atomic nucleus is fixed. Trans- 
port of X- and gamma rays takes into account Comp- 
ton scattering using the free-electron approximation, 
gamma conversion into electron-positron pair, and 
photoelectric effect. Bremsstrahlung and ionization 
are the available processes for electrons and positrons. 
The binding energy of atomic electrons is taken into 
account only at the photoelectric effect simulation. 
The simplified computational algorithm of the nu- 
clear radiation transport through the matter pro- 
vides the most calculating efficiency for the Standard 
model in comparison with other models. However, 
coherent (Rayleigh) scattering and atomic relaxation 
processes are not included for this package. 

Atomic and shell ionization effects are included 
in the Livermore model. The lower energy threshold 
of the simulation interactions is decreased down to 
250 eV. The available physical processes are Rayleigh 



and Compton scattering, photo-electric effect, pair 
production, bremsstrahlung, and ionization. The flu- 
orescence and Auger emission in excited atoms are 
also considered. This package describes the electro- 
magnetic interactions of electrons and photons tak- 
ing into account subshell integrated cross sections for 
photoelectric effect, ionization, and electron binding 
energies for all subshells. Cross-sections of particle 
interactions with matter are calculated from evalu- 
ated data libraries from Lawrence Livermore National 
Laboratory - EPDL97 for photons [11], EEDL for 
electrons [12] and EADL for fluorescence and Auger 
effects [13]. 

The physical model which originally was devel- 
oped for Penelope MC code [9] combines analyti- 
cal approach for computation of cross-section of the 
different interactions with data from cross-section 
databases calculated in Seltzer and Berger work [14]. 
Algorithm Penelope 2008 is used in GeantA begin- 
ning with version v. 4. 9. 3. This model is applicable to 
interactions in range from about 200 eV up to 1 GeV. 
Now Penelope EM Package is still tested. 

NIST data [14, 15] with empirical correction fac- 
tors for set of elements [16] are employed for compu- 
tation of bremsstrahlung loss in the EGSnrc model 
of the electromagnetic interactions. Cross-sections of 
the photoelectric effect, pair and triplet generation 
and coherent scattering cross-sections are determined 
from EDPL [11] and XCOM [17] estimations. The 
Klein-Nishina method allowing for Doppler broaden- 
ing correction and correction for chemical bond ef- 
fects is used to calculate the incoherent (Compton) 
scattering cross-sections. Parameters of atom relax- 
ation are determined for every atomic subshell. 

Contribution of the atom relaxation processes can 
be appreciable during registration of gamma-quanta 
with energies up to about 1 MeV. They are not 
described in the GeantA Standard Electromagnetic 
Model. Thus it is possible to expect the largest dif- 
ference between simulation data received with using 
Standard EM Package and EGSnrc simulation re- 
sults. 

3. RESULTS OF SIMULATIONS 

Functions of gamma-quantum energy loss distribu- 
tion (response functions) were computed for pla- 
nar detectors based on two wide band-gap semi- 
conducting compounds: mercuric (II) iodide (Hg^) 
and thallium bromide (TIBr). Hgl2 detectors are 
investigated over a long period of time [18] while 
TIBr is a comparatively new material for nuclear 
radiation detection [19]. Initial energy of gamma- 
quantum beam was chosen in the range from 0.026 up 
to 3 MeV. Simulations were run for all three GeantA 
packages of electromagnetic interactions. Geometric 
sizes of detectors and gamma-quantum beam param- 
eters were identical for all GeantA and EGSnrc mod- 
els. 

Since the thickness of investigated detectors was 
about 1 mm the minimal run length of GeantA- 
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simulated particles (cutForElectron) and gamma- 
quanta (cutFor Gamma) in detector material was set 
equal to 0.01 mm. It corresponds to the minimal en- 
ergy 3.21 keV for gamma-quanta and 43.23 keV for 
electrons in Hgl2. The minimal energy of GeantA- 
simulated gamma-quanta in TIBr is 3.85 keV and 
48.14 keV for electrons. The similar energy thresh- 
olds PCUT = 2 keV (gamma-quanta) and ECUT = 
20 keV (electrons) were set for both materials during 
EGSnrc simulations. 

In the present work we do not consider the energy 
loss of fast electrons due to the excitation of lattice 
vibrations. Statistical moments: average energy loss, 
variance and asymmetry coefficient were computed 
for every function of distribution of energy loss. Re- 
ceived dependences of the statistical moments against 
initial energy of gamma-quantum beam allow reveal- 
ing and describing in detail systematical differences 
between results of GeantA and EGSnrc simulation. 

3.1 Mercuric (II) Iodide 

Mercuric (II) iodine (Hgl2) belongs to the group of 
A IT B VI semiconductor compounds. The functions 
of distribution of gamma-quantum energy loss were 
computed for the planar Hgl2 detector with size of 
5x5x1 mm 3 . The parallel monoenergetic gamma- 
quantum beam was uniformly scanned on the detec- 
tor area. The beam incidence angle was equal to 90°. 
At the first stage of the pesent work we used GeantA 
v.4.9.3-p01 for comparison with EGSnrc simulations. 
10 6 gamma-quantum histories were simulated for ev- 
ery energy. The statistical moments of the response 
functions were calculated on their base. 

Fig. 1 shows the dependence of average en- 
ergy losses <E> of gamma-quanta (the first mo- 
ment) against initial beam energy in the investi- 
gated Hgl2 detector. Calculated energy losses are 
practically equal for all investigated models in the 
energy range less than E 1 < 0.4 MeV. The en- 
ergy losses calculated from Standard and Liver- 
more simulation data are stably lower than results 
of EGSnrc simulations in region E 1 > 0.4 MeV. 
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Fig.l. The dependence of average energy losses of 
gamma- quanta in Hgh vs initial beam energy 



The predicted mean energy losses of gamma- 
quanta received from Standard and Livermore mod- 
els practically coincides in the range between 0.4 
and 2 MeV. But the Geant4 v.4.9.3-p01 simulation 
demonstrates abrupt increase of the mean losses for 
Livermore model relative to Standard model data 
for gamma-quantum energies above E 1 > 2 MeV 
(Fig. 1). The quick increase of the electron-positron 
pair production cross-section is a single essential fea- 
ture of the gamma-quantum interaction with Hgl2 
in energy range above E 1 > 2 MeV. Consequently 
observed difference (Fig. 1) may be concerned with 
GeantA v.4.9.3-p01 errors of the positron trajectory 
simulation. 

The repeated simulation of the Hgl2-detector re- 
sponse functions using GeantA v.4.9.4-p02 (release 
from 24-06-2011) showed that these bugs have been 
fixed by authors. Moreover, as Fig. 1 shows the 
difference between Penelope model data in version 
4.9.4-p02 and result of EGSnrc simulation consid- 
erably decreased. That is why, except where noted, 
for further study we use results received from GeantA 
v.4.9.4-p02 simulations. 

To reveal the most essential differences between 
simulation results received with different GeantA EM 
packages it is necessary the detailed analysis of the 
investigated Hgl2-detector response functions in the 
energy range of gamma-quanta where different phys- 
ical processes are dominated. 

Fig. 2 shows the computed distributions of 
gamma-quantum energy losses for initial beam en- 
ergy E 1 = 0.08 MeV. At this energy the photoelec- 
tric effect is a dominating process of the gamma- 
radiation interaction with Hgl2. The amount of 
absorbed photons (on the Fig. 2 and further these 
events correspond to energy E 7 ) for all models in- 
significantly differs (Table 1). It confirms that sim- 
ilar values of photoelectric effect cross-sections for 
Hg and I are used in all models. Atomic re- 
laxation processes are not taken into account in 
the Standard model. Therefore the photoelectric 
effect event quantity is maximal for this model. 
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Fig. 2. Energy losses of gamma-quanta with energy 
E 1 = 0.08 MeV in Hgl 2 - detector 
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Table 1 and follow-up data show that the total 
efficiency of registration of gamma-quanta of investi- 
gated detectors is almost the same for all models. 

The fine structure of escape peaks is reconstructed 
from the results of EGSnrc simulations in the most 
detail. As Fig. 2 demonstrates the used energy bin 
(1 keV) allows separating KX- and LX-series of char- 
acteristic photons for both elements (Hg and I). 

Table 1. The amount of photoelectric effect events 
and total efficiency for gamma- quanta with 



0.08 MeV 



EM model 


Photoeffect 


Efficiency, Hgl2 


Standard 


804176 


0.818 


Livermore 


784069 


0.822 


Penelope 


788043 


0.822 


EGSnrc 


765231 


0.819 



In low-loss energy region of Hgl2-detector re- 
sponse functions computed with using Livermore 
and Penelope EM models show good agreement with 
EGSnrc simulation data. According to Fig. 2 these 
three models have an almost equal distribution of the 
energy losses in the region up to 0.04. . . 0.045 MeV. 
In the energy loss region above 0.045 MeV the bet- 
ter agreement with EGSnrc data is observed for 
Livermore EM model. 

All EM models lead to the similar distribution 
of the Compton tail in the gamma-quantum en- 
ergy region where Compton scattering is prevailed 
(Fig. 3). Computations with Penelope EM model 
allow obtaining the better agreement with EGSnrc 
simulation in the Compton valley whereas Stan- 
dard and Livermore results exceed EGSnrc data. 
At the same time the amount of absorbed gamma- 
quanta for the Livermore model is equal about 
two thirds of the same as other models (Table 2). 
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Fig. 3. Energy losses of gamma-quanta with energy 
E 1 =0.6 Me V in Hgl 2 - detector 

As follows Fig. 4 and Table 3 data the Livermore 
and Standard EM models employ cross-sections 
of photoelectric interaction with Hgl2 for gamma- 
quanta with energy E 1 = 1.28 MeV which values are 
much less compared with the same in the Penelope 



EM model and EGSnrc. It is confirmed by the lack 
of photopeaks for the response functions simulated 
with using these models. Moreover the disappear- 
ing of the Compton edges on the Fig. 4 indicates 
that the Livermore and Standard models employ 
smaller probabilities of the gamma-quantum Comp- 
ton scattering angels above 135° in comparison with 
the Penelope model and EGSnrc. 

Table 2. The amount of photoelectric effect events 
and total efficiency for gamma- quanta with 



E, 



0.6 MeV 



EM model 


Photoeffect 


Efficiency, Hgl2 


Standard 


12602 


0.058 


Livermore 


7418 


0.057 


Penelope 


13488 


0.058 


EGSnrc 


11312 


0.058 



If gamma-quantum energy is enough for electron- 
positron pair production (gamma-conversion) (for ex- 
ample E 1 = 1.28 MeV, Fig. 4) a double-escape 
peak of annihilation gamma-quanta becomes appre- 
ciable in the Geant4 simulations (Table 3). On 
the Fig. 4 these events correspond to the energy 
equal to 0.259 MeV. According to EGSnrc sim- 
ulation of response function of investigated Hgl2- 



detector at energy E 1 



1.28 MeV the dou- 



ble escape peak is not yet observed (Table 3). 
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Fig. 4- Energy losses of gamma- quanta with energy 
E 1 = 1.28 MeV in Hgl 2 - detector 

Table 3. The amount of the photoelectric effect and 
double escape events for gamma- quanta with 
E 1 = 1.28 MeV (Hgh) 



EM model 


Photoeffect 


Double-escape 


Standard 


39 


234 


Livermore 


18 


209 


Penelope 


1700 


201 


EGSnrc 


1083 


37 



Above results show that there is good agree- 
ment between the response functions of the investi- 
gated Hgl2-detector received with the Penelope and 
EGSnrc model apart from double-escape peak of the 
annihilation gamma-quanta. 
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Finally, the computation of the Hgl2-detector re- 
sponse function based on the 10 7 simulation trajecto- 
ries of gamma-quanta with energy E 1 = 3 MeV shows 
a good agreement between results of the Penelope 
EM model and EGSnrc (Fig. 5). Plots correspond- 
ing to results of the Livermore and Standard models 
are located appreciably lower. The escape peaks in 
that models are absent (Table 4 and Fig. 5). This 
indicates that values of cross-section of the electron- 
positron pair production employed by the Livermore 
and Standard models are smaller as compared with 
the same of the Penelope EM model and EGSnrc. 
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Fig. 5. Energy losses of gamma- quanta with energy 
E 1 = 3 MeV in Hgl 2 - detector 

Table 4- The amount of event of the escape and 
double escape of annihilation gamma-quanta (Hgh), 



E~ 



3 MeV 



EM model 


Escape 


Double-escape 


Standard 


1 


74 


Livermore 





59 


Penelope 


65 


3366 


EGSnrc 


121 


7152 



Fig. 5 shows that the response functions of the 
Hgl2-detector received with the GeantA v4.9.4-p02 
Livermore and Standard models practically coin- 
cide. So average energy losses of gamma-quanta with 
energy 3 MeV also coincide (Fig. 1). It confirms 
our supposition about the incorrect computation of 
the Hgl2-detector response functions when GeantA 
v.4.9.3-p01 Livermore and Standard models have 
been used (Fig. 1). 

Other statistical characteristics of response func- 
tions of the investigated Hgl2-detector - variance 
(Fig. 6) and skewness coefficient (Fig. 7) - also 
demonstrate a good agreement between results of 
GeantA v4.9.4-p02 Penelope EM package simulations 
and EGSnrc data. 

From Fig. 1, Fig. 6 and Fig. 7 it follows that 
main differences of all statistical characteristics of 
the Hgl2-detector response functions between GeantA 
v4.9.3-p01 and v4.9.4-p02 simulation data correspond 
to the electron-positron pair production region (E 1 



> 1.022 MeV). In the gamma-quantum energy re- 
gion where positrons are not generated the response 
functions received with different Geant4 versions 
coincide. Therefore the cause of anomalous simu- 
lation results is incorrect values of specific ioniza- 
tion energy losses of positrons in Hgl2 those used 
by GeantA v4.9.3-p01. The differences between re- 
sponse functions of the investigated TIBr-detector re- 
ceived using of both GeantA versions are insignificant. 
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Fig. 6. The dependence of variance of the Hg^- 
detector response functions vs gamma- quantum 
energy 
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Fig. 7. The dependence of skewness coefficient of the 
Hgl2-detector response functions vs gamma-quantum 
energy 

3.2 Thallium Bromide 

Thallium bromide (TIBr) belongs to the group of 
pjii^QVi semiconductor compounds. Quality differ- 
ences between the linear gamma-ray attenuation co- 
efficients of TIBr and Hgl2 are negligible. These dif- 
ferences are characteristic for the energy region where 
photoelectric effect is dominated (Fig. 8). At the 
same time the TIBr linear gamma-ray attenuation 
coefficient exceeds the same of Hgl2 about 10% be- 
ginning from gamma-ray energy E 1 =0.1 MeV. As 
follows we can expect that average losses of gamma- 
quanta energy in TIBr will be higher than in Hgl2 at 
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the same detector thickness. 

Distributions of the gamma-quantum energy 
losses were calculated for the planar TIBr-detector 
with size of 3.142 mm 2 x 0.8 mm. The parallel mo- 
noenergetic beam of gamma-quanta was uniformly 
distributed on the detector area. The beam incidence 
angle was equal to 90°. 10 6 histories were simulated 
for each gamma-ray energy. Statistical moments of 
response functions were calculated using these histo- 
ries (Fig. 9 and Fig. 10). 

The data plotted on the Fig. 9 - Fig. 12 are re- 
ceived from the GeantA v4.9.4-p02 simulations. Re- 
sponse functions of the investigated TIBr-detector 
obtained using version v. 4. 9. 3 coincide with ver- 
sion v.4.9.4 calculations within the scope of statis- 
tical uncertainties. It gives grounds for affirma- 
tion that the reason of above-mentioned anoma- 
lous simulation results for Hgl2-detector is not the 
errors of the positron trajectory simulation algo- 
rithm. As we think these anomalous results are con- 
sequence of incorrect values of the specific ioniza- 
tion energy losses of positrons in Hgl2. Similarly 
the case of Hgl2-detector the better coincidence be- 
tween EGSnrc and GeantA simulation data is ob- 
served when the Penelope model of electromagnetic 
interactions has been used (Fig. 9 and Fig. 10). 




E , MeV 

Fig. 8. The linear gamma-ray attenuation coefficient 
for Hgl 2 and TIBr [17] 

The analysis of simulated TIBr-detector response 
functions for gamma-ray photons with energy E 1 = 
0.6 MeV (Fig. 11) shows the presence of escape peaks 
of characteristic photons of Tl and Br KX-series in 
all models apart from the Standard EM model. The 
Compton tail is identically reconstructed in all sim- 
ulations. As for the Compton valley region that the 
better agreement with EGSnrc data is obtained for 
the GeantA Penelope EM package simulation. The 
amount of the absorbed gamma-quanta with the en- 
ergy 0.6 MeV for the Livermore EM model simu- 
lation is about 30% less in comparison with other 
models (Table 5). 

The response function of the investigated TIBr- 
detector for rather high energy gamma-quanta (E 1 
= 3 MeV, Fig. 12) received with using EGSnrc 
MC package shows a good agreement only with 



the results of the GeantA Penelope model simu- 
lation. For the Livermore and Standard models 
the escape peak of the annihilation gamma-quanta 
is almost absent and amplitude of the double es- 
cape peak is more than an order of magnitude 
less as compared with other simulations (Table 6). 




E , MeV 

Fig. 9. The dependence of average energy losses of 
gamma- quanta in TIBr vs the initial energy beam 
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Fig. 10. The dependence of variance of TIBr- 
detector response functions vs the gamma- quantum 
energy 
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E, MeV 

Fig. 11. The energy losses of gamma-quanta with 
the energy E 1 =0.6 MeV in TIBr-detector 
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Table 5. The amount of photoelectric effect events 
and total efficiency for gamma- quanta with 



E 1 = 0.6 MeV 



EM model 


PhotoefTect 


Efficiency, TIBr 


Standard 


14288 


0.061 


Livermore 


9231 


0.060 


Penelope 


15652 


0.060 


EGSnrc 


16804 


0.064 
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Fig. 12. The energy losses of gamma- quanta with 5. 
the energy E 1 = 3 MeV in TIBr- detector 

Table 6. The amount of event of the escape and 
double escape of annihilation gamma-quanta (TIBr), 

E 1 = 3 MeV 6 ' 



EM model 


Escape 


Double-escape 


Standard 


13 


435 


Livermore 


5 


395 


Penelope 


235 


5881 


EGSnrc 


577 


13433 



4. CONCLUSIONS 

The comparison of the response functions of the Hgl2 
and TIBr gamma-radiation detectors computed for 
all models of the electromagnetic interactions of the 
GeantA package and EGSnrc has been made in this 
work. In the energy region of gamma-quanta up 
to 0.4 MeV there is a good agreement between re- 
sponse functions of the investigated detectors cal- 
culated with using EGSnrc and both Livermore 
and Penelope models of the GeantA package. For 
wider gamma-quantum energy region between 0.026 
and 3 MeV statistical parameters of the response 
functions simulated with EGSnrc and GeantA pack- 
ages have a good agreement only when the Penelope 
model of electromagnetic interactions has been used. 

By the Hgl2-detector example it is shown that the 
comparison of the response functions received with 
using different simulation packages allow removing 
incorrect results and explaining the reason of their 
appearance. 
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